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Abstract 

The fundamental problem of non-singular dislocations in the framework of the 
theory of gradient elasticity is presented in this work. Gradient elasticity of Helmholtz 
type and bi-Helmholtz type are used. A general theory of non-singular dislocations 
is developed for linearly elastic, infinitely extended, homogeneous, and isotropic 
media. Dislocation loops and straight dislocations are investigated. Using the the- 
ory of gradient elasticity, the non-singular fields which are produced by arbitrary 
dislocation loops are given. 'Modified' Mura, Peach-Koehler, and Burgers formu- 
lae are presented in the framework of gradient elasticity theory. These formulae 
are given in terms of an elementary function, which regularizes the classical ex- 
pressions, obtained from the Green tensor of the Helmholtz-Navier equation and 
bi-Helmholtz-Navier equation. Using the mathematical method of Green's func- 
tions and the Fourier transform, exact, analytical, and non-singular solutions were 
found. The obtained dislocation fields are non-singular due to the regularization of 
the classical singular fields. 
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1 Introduction 

Dislocations play an important role in the understanding of many phenomena in solid 
state physics, materials science, and engineering. They are the primary carriers of crystal 
plasticity. Dislocations are line defects which can be straight or curved lines. The internal 
geometry of generally curved dislocations, in deformed crystals is very complex. In the 
classical theory of dislocation loops in isotropic materials (DeWit, 1960; Lardner, 1974; 
Hirth and Lothe, 1982) two key equations are the Burgers formula (Burgers, 1939) for 
the displacement, and the Peach-Koehler formula (Peach and Koehler, 1950) for the 
stress. These equations are very important for the interaction between complex arrays 
of dislocations. The classical description of the elastic fields produced by dislocations is 
based on the theory of linear elasticity. There is the problem of mathematical singularities 
at the dislocation core, and an arbitrary core-cutoff radius which must be introduced to 
avoid divergence. In the classical continuum theory of dislocations (Kroner, 1958; DeWit, 
1960; Nabarro, 1967; Lardner, 1974; Mura, 1987; Teodosiu, 1982; Li and Wang, 2008) the 
concept of Volterra dislocations is used, and the dislocation core is described by a Dirac 
delta function. This unsatisfactory situation can only be remedied, when the fact that 
physical dislocations have a finite core region and no singularities exist are taken into 
account. 

As already pointed out by Kroner (1958) and Lothe (1992) the divergence can be 
avoided when dislocation distributions other than the delta function are used. Lothe 
(1992) considered a standard core model with constant density of dislocations in a planar 
strip with width d. However, the value of d remains undetermined and the expressions 
for the elastic fields are more complicated than their singular counterparts and difficult to 
use for generally curved dislocations. Moreover, for non-planar configurations the theory 
becomes much more complex. 

In order to remove the singularities of dislocations and to model the dislocation core 
more realistically, continuum theories of generalized elasticity may be used. A very promis- 
ing candidate of such a theory is the so-called gradient elasticity theory. The theory of 
gradient elasticity was originally proposed by Mindlin (1964, 1965), and Mindlin and 
Eshel (1968) (see also Eshel and Rosenfeld (1970)). The correspondence between the 
strain gradient theory and the atomic structure of materials with the nearest and next 
nearest interatomic interactions was exhibited by Toupin and Grazis (1964). The origi- 
nal Mindlin theory possesses too many new material parameters. For isotropic materials, 
Mindlin's theory of first strain gradient elasticity (Mindlin, 1964; Mindlin and Eshel, 1968) 
involves two characteristic lengths, and Mindlin's theory of second strain gradient elastic- 
ity (Mindlin, 1965) possesses four characteristic lengths. The discrete nature of materials 
is inherently incorporated in the formulations through the characteristic lengths. The 
capability of strain gradient theories in capturing size effects is a direct manifestation of 
the involvement of characteristic lengths. Simplified versions, which are particular cases 
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of Mindlin's theories, were proposed and used for dislocation modelling. Such simplified 
gradient elasticity theories are known as gradient elasticity of Helmholtz type (Lazar and 
Maugin, 2005), with only one material length scale parameter and gradient elasticity of 
bi-Helmholtz type (Lazar and Maugin, 2006a; Lazar et al., 2006a) which involves two ma- 
terial length scale parameters as new material coefficients. Gradient elasticity is a contin- 
uum model of dislocations with core spreading. Non-singular fields of straight dislocations 
were obtained in the framework of gradient elasticity of Helmholtz type by Gutkin and 
Aifantis (1999); Lazar and Maugin (2005, 2006a); Lazar et al. (2005) and Gutkin (2000, 
2006) (see also, Gutkin and Ovid'ko (2004)). Surprisingly enough up until now, not a 
single work has been done in the direction of non-singular dislocation loops using strain 
gradient elasticity theory. The reason may be in the expected mathematical complexity 
of the problem. Such non-singular solutions of arbitrary dislocation loops could be very 
useful for the so-called discrete dislocation dynamics (e.g., Li and Wang (2008); Ghoniem 
and Sun (1999)). 

The aim of this paper is to present non-singular solutions of arbitrary dislocation 
loops, by using simplified gradient elasticity theories. We present the key-formulae of 
dislocations loops valid in the framework of gradient elasticity, and also reemphasize 
straight dislocations in gradient elasticity. The technique of Green functions for the key- 
formulae is used, and analytical closed-form solutions for the dislocation fields are derived. 

The paper is organized as follows. In Section 2, the fundamentals of gradient elasticity 
of Helmholtz type are given. Dislocation loops and straight dislocations are investigated. 
In Section 3, the theory of gradient elasticity of bi-Helmholtz type is considered. Dislo- 
cation loops and straight dislocations will be examined in this framework. In Section 4, 
the conclusions are given. All the mathematical and technical details are given in the 
Appendices. 

2 Gradient elasticity of Helmholtz type 

A straightforward framework to obtain non-singular fields of dislocations is the so-called 
theory of gradient elasticity. A simplified theory of strain gradient elasticity is called gradi- 
ent elasticity of Helmholtz type (Lazar and Maugin, 2005, 2006a). This gradient elasticity 
of Helmholtz type is a particular gradient elasticity theory evolving from Mindlin's general 
gradient elasticity theory (Mindlin, 1964; Mindlin and Eshel, 1968). This theory is also 
known as dipolar gradient elasticity theory (Georgiadis, 2003), simplified strain gradient 
elasticity theory (Gao and Ma, 2010a,b) and special gradient elasticity theory (Altan and 
Aifantis, 1997). The theory of gradient elasticity of Helmholtz type is the gradient version 
of Eringen's theory of nonlocal elasticity of Helmholtz type (Eringen, 1983, 2002) which 
is well-established. 

The strain energy density of such a simplified gradient elasticity theory for an isotropic, 
linearly elastic material has the form (Lazar and Maugin, 2005; Gao and Ma, 2010a) 

W — - Cijki(3ij(3ki + - £ 2 Cijkid m (3ijd m (3ki , (1) 
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where Cijki is the tensor of elastic moduli with the symmetry properties 

Cijki = Ckuj = Cjiki = Cijik (2) 
and it reads for an isotropic material 

C^ki = ^(Sik^ji + SuSjk) + A SijSki , (3) 

where fj, and A are the Lame moduli, denotes the elastic distortion tensor. If the 
elastic distortion tensor is incompatible, it can be decomposed as follows 

Pij = d jUi - Pl , (4) 

where Ui and (3fj denote the displacement vector and the plastic distortion tensor, re- 
spectively. In addition, i is the material length scale parameter of gradient elasticity 
of Helmholtz type. For dislocations, I is related to the dislocation core radius and is 
proportional to a lattice parameter. Due to the symmetry of C^ki, Eq. (1) is equivalent 
to 

X x 

W = - Cij k ieijeki + - £ Cij k [d m eijd m e k i , (5) 

where = 1/2 + is the elastic strain tensor. The condition for non- negative 
strain energy density, W > 0, gives 

(2/i + 3A)>0, /i>0, £ 2 >0. (6) 

The reason that the elastic and plastic distortion tensors are incompatible can be 
the presence of dislocations. Dislocations cause self-stresses that means stresses caused 
without the presence of body forces. The dislocation density tensor is defined in terms of 
the elastic and plastic distortion tensors as follows (e.g., Kroner (1958)) 

otij = ejkidkPu (7) 
«»j = -£jkid k (3fi (8) 

and it fulfills the Bianchi identity of dislocations 

djOij = , (9) 

which means that dislocations do not end inside the body. Eq. (9) is a 'conservation' law 
and shows that dislocations are source-free fields. 

From Eq. (1) it follows that the constitutive equations are 

dW dW 

&ij — -QQ- — 7^— = Cijkiftki = CijkiCki , (10) 
dW dW 2 2 

Tijk oo n ^ Cijmndkfimn £ &k®ij i (H) 

COkPij OOkCij 



were are the components of the Cauchy stress tensor, are the components of the 
so-called double stress tensor. It can be seen that t is the characteristic length scale for 
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double stresses. Using Eqs. (10) and (11), Eq. (5) can also be written as (Lazar and 
Maugin, 2005) 

W = ^<T ij e ij + -£ 2 d k a ij d k e ij . (12) 

The strain energy density (12) exhibits the symmetry both in and and in dk&ij and 

The total stress tensor is given as a combination of the Cauchy stress tensor and the 
divergence of the double stress tensor 

a ij = a ij - dkTijk = (1 - (13) 

and it fulfills the equilibrium condition for vanishing body forces 

d j°ij = d A a ij ~ d kn jk ) = . (14) 

The stress tensor cr?- is called in the notation of Jaunzemis (1967) the polarization of the 
Cauchy stress a^. Due to gradient elasticity of Helmholtz type, Eq. (13) reduces to an 
inhomogeneous Helmholtz equation where the total stress tensor is the inhomogeneous 
piece. As pointed out by Lazar and Maugin (2005, 2006a), the total stress tensor may be 
identified with the singular classical stress tensor. This identifies that the inhomogeneous 
Helmholtz equation (13) is in full agreement with the equation for the stress proposed 
by Eringen (1983, 2002) in his theory of nonlocal elasticity of Helmholtz type. 

As shown by Lazar and Maugin (2005, 2006a) the following governing equations for 
the displacement vector, the elastic distortion tensor, the dislocation density tensor, and 
the plastic distortion tensor can be derived in the framework of gradient elasticity of 
Helmholtz type 



where 



Lu t = u°, (15) 

= (16) 

La ij = a° ij , (17) 

Lf% = f%°, (18) 



L = l-fA (19) 

is the Helmholtz operator. The singular fields u®, and are the sources of the 

inhomogeneous Helmholtz equations (15)-(18). The Helmholtz equations (15) and (16) 
can be further reduced to Helmholtz-Navier equations 

L LikUk = Cijkidj(3^i° , (20) 

L Lihfikm — —Cijkltmlrdja^ , (21) 

where Lik = Cijkidjdi is the differential operator of the Navier equation. For an isotropic 
material, it reads 

L ik = fJ> 5ifeA + (fj, + X) did k . (22) 
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In Eqs. (20) and (21) the sources are now the plastic distortion (3^° and the dislocation 
density a° kr known from classical elasticity. 

The corresponding three-dimensional Green tensor of the Helmholtz-Navier equation 
is defined by 



L L ik G kj = —Sij S(x-x') 
and is calculated as (see Eq. (B.14)) 



Gij(R) 



167TyU(l — v) 



with 



2(1 - u)SijA - didj A(R) , 



A(R) = R+^-(l-e-^ e ) 



(23) 



(24) 



(25) 



and R — \x — x'\. In the limit I — > 0, the three-dimensional Green tensor of classical 
elasticity (Mura, 1987; Li and Wang, 2008) is recovered from Eqs. (24) and (25). It is 
important to note that A(R) can be written as the convolution of R and G(R): 



A(R) = R* G(R) 



(26) 



where * denotes the spatial convolution and G is the three-dimensional Green function of 
the Helmholtz equation 



1 

AttPR ' 



LG = 5(x - x') . 

It reads (Wladimirow, 1971) 

G(R) = 

In addition, it holds 

AAR = 

The function (25) fulfills the relations 

LAAA(R) 



R/e 



-8tt8(x — x') . 



-87T 5(x — x') , 
AAA(R) = -8irG(R), 
LA(R) = R. 



(27) 



(28) 



(29) 



(30) 
(31) 
(32) 



Thus, A(R) is the Green function of Eq. (30) which is a Helmholtz-bi-Laplace equation. 

Using Eqs. (A. 3) and (A. 4) for the differentiation of Eq. (24), the explicit form of the 
three-dimensional Green tensor of the Helmholtz-Navier equation is obtained 



Gij{R) 



1 



167TyU(l — V) 



^ ((3 - Av) (l - e- R /^ + -1 (2£ 2 - (R 2 + 2£R + 2£ 2 )e- R ^ 



+ 



RiRj 



1 - 



6^ 
R? 



+ 2 + 



R ' R 2 )' 



6£ 6£ : 
+ 



-n/e 



(33) 
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which is non-singular. It is worth noting as a check, that Eq. (33) is in agreement with 
the corresponding expressions derived by Polyzos et al. (2003) and Gao and Ma (2009) 
using slightly different approaches. 

Note, that the Green tensor (33) gives the non-singular displacement field, U{ = Gijfj 
(fj is the constant value of the magnitude of the point force acting at the arbitrary 
position x' in an infinite body), of the Kelvin point force problem (e.g., Gurtin (1972); 
Mura (1987); Hetnarski and Ignaczak (2004)) in the framework of gradient elasticity of 
Helmholtz type. The original solution of a concentrated force in an infinite body in the 
context of the classical continuum theory of elasticity was given by Lord Kelvin (1882). 

2.1 Dislocation loops 

In this subsection, the characteristic fields of dislocation loops in the framework of gradient 
elasticity theory of Helmholtz type are calculated. 

For a general (non-planar or planar) dislocation loop L, the classical dislocation density 
and the plastic distortion tensors are (e.g., DeWit (1973a); Kossecka (1974)) 

a% = b t 5j (L) =bij> 5(x- x') dL' 3 , (34) 

P%° = ~ b i S i(S) = ~bi I S{x - x') dS'j , (35) 

Js 

where hi is the Burgers vector of the dislocation line element dL'j at x' and dSj is the 
dislocation loop area. The surface S is the dislocation surface, which is a cap of the 
dislocation line L. Sj(L) is the Dirac delta function for a closed curve L and Sj(S) is the 
Dirac delta function for a surface S with boundary L. 

The solution of Eq. (17) can be written as the following convolution integral 

oiij = G * a% = hj G(R) dL; , (36) 

where G(R) denotes the three-dimensional Green function of the Helmholtz equation 
given by Eq. (28). The explicit solution of the dislocation density tensor for a dislocation 
loop in gradient elasticity is calculated as 

^iM^r^" (37) 

describing a spreading dislocation core distribution. The plastic distortion tensor of a 
dislocation loop, which is the solution of Eq. (18), is given by the convolution integral 

$ij — G* /f° = -h [ G(R) dS'j . (38) 
Js 

It reads as 

h. r e~ R / £ 

(39) 
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Substituting Eq. (39) in Eq. (8) and using the Stokes theorem, we obtain formula (37). 

Using the Green tensor (24), and after a straightforward calculation all the generaliza- 
tions of the Mura, Peach-Koehler, and Burgers formulae towards gradient elasticity can 
be obtained. Starting with the elastic distortion tensor of a dislocation loop, the solution 
of Eq. (21) gives the representation as the following convolution integral 



e mnr CjkinGij,k(R)o^ r (x') dV' , 



(40) 



where Gij^ = dkGij. Substituting the classical dislocation density tensor of a dislocation 
loop (34) and carrying out the integration of the delta function, we find the modified 
Mura formula valid in gradient elasticity 



A, 



i( X ) — j> e mnrblCjkl n Gij t k(R) dL' r . 



(41) 



Substitute Eqs. (3) and (24) into Eq. (41) and obtain after rearranging terms 



8tt 



(hd n - brA + biS in di)A + j— (b n diA - hdidnd^ A{R) dL' r . 



(42) 



Using the identity 



[trstfiij — £jst8ir)b s dt 



(43) 



€rjn(bid n b n di) ^rjn^kin^kstbsdt ($rk$ji $ri$ jkj^kstbsdt 
= {trkl$ij — £jkAr)bkdl 

and the relation 

<f> e rjn {b n di - bid n )didiA(R) dL' r = <h b s (e jst d r - e r8t dj)d t d i A{I{) dL' r 

JL JL 

= <b d(b s e jst d t diA(R)) - (t b s e rst d t d j d i A(R)dL' r = - <p 6 s e rat ^9 7 -9 i >l(i2) dL' r 
Jl Jl Jl 

= -£b k e rkl d l d j d i A(R)dL' r , (44) 

the non-singular elastic distortion (42) of a dislocation loop becomes 

Pij(x) = j> (e jk iS ir - e rk i5 i:j + e rij 5 M )diA + — — e^d^d,] A(R) dL' r . (45) 

This is the 'Mura formula' for a dislocation loop in gradient elasticity. It is important to 
note that if Eq. (45) is substituted into (7) and the relation (31) is used, the dislocation 
density of a dislocation loop (37) is recovered. 

The symmetric part of the elastic distortion tensor (45) gives the elastic strain tensor 
of a dislocation loop 

eij(x) = j> e jk iS ir + ^ e ik i5 jr - e rfc z<%)<9zA + e rkX dididj A(R) dL' r . (46) 
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The elastic dilatation of a dislocation loop is nothing but the trace of the elastic strain (46) 

(1 - 2u)b k 



8tt( 



I ^ jf e rkl d^A{R) dL' r 



(47) 



The elastic rotation vector is defined as the skewsymmetric part of the elastic distortion 
tensor oj\ = ^tijiPij and reads 



u^x) = ~^j (hrdk - \ 6 kr di)AA(R) dL' r 



(48) 



Using the constitutive relation (10) with Eq. (3), the non-singular stress field produced 
by a dislocation loop is found 

an (x) = j> (e jkl 5 ir + e m 5 jr ) d t A + j— - e rkl (didj - <J - A) d t A(R) dL' r , (49) 

which can be interpreted as the Peach-Koehler formula within the framework of gradient 
elasticity. One may verify that the stress is divergence-less, djO^ = 0. The double stress 
tensor of a dislocation loop is easily obtained if Eq. (49) is substituted into Eq. (11). 
The solution of Eq. (20) is the following convolution integral 



/oo 
CjkmGij^Pffix') dV. 
-oo 



(50) 



Substituting the classical plastic distortion of a dislocation loop (35) into Eq. (50) gives 
the modified Volterra formula valid in gradient elasticity 



Ui(x) = / biC jtan G ijJe (R)dS' n . (51) 

Substituting Eqs. (3) and (24) into Eq. (51) and rearranging terms yield 

b f I 

Ui (x) = / / (6ud n + 5 in dt - 5i n di)A + (5 ln A - did n )di A(R) dS' n . (52) 



Except the first term of Eq. (52), we apply the Stokes theorem in order to obtain line 
integrals with 



and 



[ M - 8i n di)AA(R) dS' n = -<f e rtl AA(R) dL' r 
Js Jl 

(S nl A - did^diAiR) dS' n = -j erijdjdiA{R) dL' r . 



(53) 



(54) 



In this way, the key-formula for the non-singular displacement vector in gradient elasticity 
is found 



Ui{x) 



8tt 



AdjA(R) dS' + 



bl^rlj 

~8^T 



SijA - 



1 - v 



did, \ A(R) dL' r , (55) 
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which is the Burgers formula in the framework of gradient elasticity of Helmholtz type. 
Eq. (55) determines the displacement field of a single dislocation loop. The Eqs. (45)- 
(55) are straightforward, simple, and closely resemble the singular solutions of classical 
elasticity theory. In the limit i — >■ 0, the classical expressions are recovered in Eqs. (45)- 
(55). The expressions (45), (49), and (55) retain most of the analytic structure of the 
classical Mura, Peach-Koehler, and Burgers formulae. The expressions (45)-(55) are 
given in terms of the elementary function A(R) given in Eq. (25), instead of the classical 
expression R. The explicit expressions can be obtained by simple substitution of the 
formulae for the derivatives of A given in Eqs. (A.2)-(A.6). It is important to note 
that Eqs. (45)-(55) are non-singular due to the regularization of the classical singular 
expressions (see Appendix A). As an example, we substitute Eqs. (A. 5) and (A. 6) into 
Eq. (49) and obtain the explicit expression for the stress tensor 



2 



£rkl< 



\2Rt 

7 i?3 



1 + 



jY 



-R/e 



+ 



Sij Rk + Rj + $jk Ri 



1 - V 

3RiRjRk 



R 3 



6f 
R 2 



-R/e 



r 2-| 



-R/e 



10f 



i 



-R/e 



A m 2R\ 

4 + — + — 1 e 



R 



3£J 



R/e 



dL' r . (56) 



To give the expression (55) more explicitly. Using Eq. (A. 6) we introduce a generalized 
solid angle valid in gradient elasticity of Helmholtz type 

il(x,t) = -\ f^d 3 A{R)<lS' 3 = / 5 §(l-(l + f)e-^)d^. (57) 

Eq. (57) is non-singular and depends on the length scale £. In the limit £ — > 0, the usual 
solid angle (e.g., Li and Wang (2008)) is recovered. Thus, using Eq. (57) and carrying 
out some differentiations with Eqs. (A. 2) and (A. 4), we obtain from Eq. (55) the explicit 
gradient elasticity version of the Burgers formula 



uAx) 



An 



8tt(1 - v) 



The simplicity of our results is based on the use of gradient elasticity theory of Helmholtz 
type. Our results can be used in computer simulations of dislocation cores at nano-scale 
and in numerics as fast numerical sums of the relevant elastic fields as it is used for the 
classical equations (e.g., Ghoniem and Sun (1999)). 



2.2 Straight dislocations 

In this subsection, using the modified Mura equation of gradient elasticity of Helmholtz 
type (40), the non-singular elastic distortion fields of straight dislocations as a check of 
our general approach are calculated. 
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2.2.1 Screw dislocation 



A screw dislocation corresponds to the anti-plane strain problem. The Green function 
of the anti-plane strain problem in gradient elasticity of Helmholtz type is nothing but 
the Green function of the two-dimensional Helmholtz-Laplace equation and it reads (see 
Eq. (B.21)) 

G ZZ {R) = ~^{lE + lnR + K (R/£) } , (59) 

where R = yjjx — x') 2 + (y — y') 2 , ^ E is the Euler constant and K n is the modified Bessel 
function of order n. The Green function (59) is non-singular. The gradient of the Green 
function (59) is obtained as 



G 



Next, substituting Eq. (60) and the dislocation density of a screw dislocation a Q zz = 
b z 5(x)8(y) into Eq. (40), the elastic distortion produced by a screw dislocation is obtained. 
For an infinite screw dislocation along the z-axis with Burgers vector b z , the non-singular 
components for the elastic distortion are calculated as 



b z y 



2tt r 
b z x r r 



[l-jKMt)}, (61) 



zy 2tt r 



where r = \J x 2 + y 2 . The expressions obtained earlier by Lazar (2003) and Lazar and 
Maugin (2006a) are recovered. In the limit £ — > 0, the classical expressions given by DeWit 
(1973b) are recovered in Eqs. (61) and (62). 

The Green function (59) gives the non-singular displacement field u z = —G zz f z of a 
line force with the magnitude f z calculated by Lazar and Maugin (2006b) in the framework 
of gradient elasticity. 

2.2.2 Edge dislocation 

Now the plane strain problem of an edge dislocation is investigated. The Green tensor 
of the plane strain problem in gradient elasticity of Helmholtz type is derived as (see 
Eq. (B.17)) 



Gij(R) = Sij^E + lnR + K (R/£) } 



1 

4 



167r//(l — v) 



-—- d t dj{R 2 { lE + In R) + U 2 { lE + \nR + K (R/£) ) } . (63) 



It is obvious that the terms proportional to the Euler constant do not contribute to the 
elastic distortion fields. The two-dimensional Green tensor (63) is non-singular. In the 
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limit £ — > 0, the two-dimensional Green tensor of classical elasticity (Mura, 1987; Li and 
Wang, 2008) is recovered in Eq. (63). The gradient of the Green tensor (63) is given by 



GijA R ) 



8tt/z(1 - v) 



to A \ X ^ k X X , O RiRjRk 

(3-4,)^-^-^^-^^ + 2^^ 



+ ^ { 5 ^ Rk + **** + ^ - 4 ^f^) - K " W) ) 



(64) 



Substituting Eq. (64) and the dislocation density of an edge dislocation along z axis 
with Burgers vector b x , a° xz = b x S(x)S(y), into Eq. (40), the elastic distortion of an edge 
dislocation is obtained. Eventually, the non- vanishing components of the elastic distortion 
of an edge dislocation are calculated as 



b x V f,, n ^ 2x 4 ^ /2 o 2 x 2(y 2 -3x 2 ) 

_2jS-^ Ki(rm y (65) 



{<3-*>-£-^-V) + ^^(r/<) 

A'tfr/C)), (66) 



47r(l — z/) r 

2 (y2 + (l- y ) r 2) 

x r, t , , 2y 2 , A£\ 2 2 . 2(x 2 -3y 2 ) 



" 4tt(1 - !/) r 2 



+ 2 (, 2 -(i-^) A , (r/o}| (g7) 



= ~ 4tt(i-i/) ^ l ( " } " T 2 " " 7^ (y " } + 7* 2(r/ } 



2 (a; 2 — vr 
Tr 



Ki(r/£)}, (68) 



which are non-singular and agree with the formulae given by Lazar (2003) and Lazar and 
Maugin (2006a). In the limit I — > 0, we obtain in Eqs. (65)-(68) the classical expressions 
given by DeWit (1973b). As discussed by Lazar et al. (2006a), the dislocation core radius 
can be defined straightforwardly in the framework of gradient elasticity as R c ~ 6£. If 
£ ~ 0.4 a, where a denotes the lattice parameter, is adopted as proposed by Eringen 
(1983), the dislocation core radius is R c ~ 2.5a. Using £ ~ 0.4a, the internal length 
reduces to £ ~ 1.97 A for lead (Pb) with a = 4.95 A. 

Note that the two-dimensional Green function (63) gives the non-singular displacement 
field, Ui = —Gijfj, of a line force with magnitude fj calculated by Lazar and Maugin 
(2006b) in the framework of gradient elasticity. 
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3 Gradient elasticity of bi-Helmholtz type 



In this section, gradient elasticity theory of higher order is considered. Gradient elasticity 
theory of higher order was originally introduced by Mindlin (1965, 1972) (see also, Jaun- 
zemis (1967); Wu (1992); Agiasofitou and Lazar (2009)). Mindlin's theory of second strain 
gradient elasticity involves for isotropic materials, in addition to the two Lame constants, 
sixteen additional material constants. These constants produce four characteristic length 
scales. 

A simple and robust gradient elasticity of higher order which is called gradient elas- 
ticity theory of bi-Helmholtz type was introduced by Lazar et al. (2006a) and Lazar and 
Maugin (2006a) and successfully applied to the problems of straight dislocations (Lazar 
et al., 2006a; Lazar and Maugin, 2006a), straight disclinations (Deng et al., 2007) and 
point defects (Zhang et al., 2006). Lazar et al. (2006a) and Lazar and Maugin (2006a) 
have shown that all state quantities are non-singular. By means of this second order 
gradient theory it is possible to eliminate not only the singularities of the strain and 
stress tensors, but also the singularities of the double and triple stress tensors and of the 
dislocation density tensors of straight dislocations at the dislocation line. In general, all 
fields calculated in the theory of gradient elasticity of bi-Helmholtz type are smoother 
than those calculated by gradient elasticity theory of Helmholtz type. In general, there a 
two main motivations for the use of gradient elasticity of bi-Helmholtz type: a consistent 
regularization of all state quantities, and a more realistic modelling of dispersion relations. 
A simple higher-order gradient theory in order to investigate dislocation loops should be 
used. The theory of gradient elasticity of bi-Helmholtz type is the gradient version of 
nonlocal elasticity of bi-Helmholtz type (Lazar et al., 2006b). 

The strain energy density of gradient elasticity theory of bi-Helmholtz type for an 
isotropic, linearly elastic material has the form (Lazar et al., 2006a) 

1 1 2 1 4 

W — - Cijki(3ij(3ki + - £iCijkid m (3ijd m (3ki + - £ 2 Cijkidnd m (3ijd n d m (3ki , (69) 

where l\ = £, £ 2 is another characteristic length scale and C^i is given in (3). Due to the 
symmetry of Cijki, Eq. (69) is equivalent to 

W = — C \jkl^ijCkl + — ^iCijkldm e ijd m e kl + ^ ^2^ \jkld n dmZijd n dm e kl ■ (70) 

In addition to the constitutive equations (10) and (11) another one is present in such a 
higher-order gradient theory, 

dW dW 4 4 

r ijkl = oo o n = o = ^2 Cij mn dldk(3 mn = £ 2 9;9fcC7jj , (71) 

OOiOkPij OOiOk^ij 

where t^h is called the triple stress tensor. It can be seen that £2 is the characteristic 
length scale for triple stresses. On the other hand, l\ is the characteristic length scale for 
double stresses. Using Eqs. (10), (11), and (71), Eq. (70) can also be written as (Lazar 
et al, 2006a) 

11 1 

w = 2 Gi 3 e v + 2 ^(Tijdkeij + - ijdidkaijdidkeij . (72) 
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The strain energy density (72) exhibits the symmetry in and e^, in d k a %j and d k eij, 
and in d\d k Oij and did^e^. The condition for non-negative strain energy density, W > 0, 
gives 

^ > , 4 > , (73) 

in addition to (3/i + 2A) > and /x > 0. 
The total stress tensor reads now 

a ij = a ij ~ dkTijk + didkTijkl ■ 

In absence of body forces, the equation of equilibrium has the following form 

dj a ij = dj( a ij ~~ 9 k T ijk + did k T ijk i) = . 
Using Eqs. (11) and (71), the total stress tensor (74) can be written 

where the differential operator L is given by 

L = (1 - £jA + £ 4 2 AA) = (1 - c\A) (1 - cjA) 

with 
and 

£l = cl + cl, (80) 
t\ = c\c\. (81) 

Due to its structure as a product of two Helmholtz operators, the differential operator (77) 
is called bi-Helmholtz operator. 

An important point, is the question concerning the mathematical character of the 
length scales c\ and c 2 . Mindlin (1965) (see also, Mindlin (1972); Wu (1992)) pointed 
out that the conditions for non-negative W supply no indications of the character, real 
or complex, of the characteristic lengths. Mindlin (1965) and Wu (1992) have treated 
the characteristic lengths as if they were real and positive. They also pointed out that a 
complex character of the lengths is equally admissible. The character, real or complex, 
of the lengths dictates the behaviour of the field variables. In the theory of gradient 
elasticity of bi-Helmholtz type the condition for the character, real or complex, of the 
length scales c\ and c 2 can be obtained from the condition if the argument of the square 
root in Eqs. (78) and (79) is positive or negative. Thus, c\ and c 2 are real if 

t\ - U\ > , (82) 
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(74) 
(75) 
(76) 
(77) 

(78) 
(79) 



and ci and c 2 are complex if 

i\ - U\ < . (83) 

If the lengths C\ and c 2 are complex, then the behaviour of the solutions of the field 
quantities would be oscillatory. In this case, the far-field behaviour of the strain and 
stress fields of dislocations would not agree with the classical behaviour. The limit from 
gradient elasticity of bi-Helmholtz type to gradient elasticity of Helmholtz type is: c 2 — > 0, 
£ 2 — )> and ci — > l\. If ci is complex, then also t\ becomes complex what would be rather 
strange. Thus, a real character of the length scales c\ and c 2 seems to be more realistic and 
more physical. In addition, Zhang et al. (2006) determined, in an atomistic calculation, 
the length scales c\ and c 2 as positive and real for graphene. In what follows, the length 
scales Ci and c 2 will be treated as if they are real and positive. 

The Green tensor of the bi-Helmholtz-Navier equation is calculated as (see Eq. (B.27)) 

G ^ R) = 167^(1 - y) t 2(1 " u)5ijA " dld ^ A(R) ' (M) 

where the elementary function (25) is changed to 

A(R) = R + - — ^ 1 (4 e-*« - 4 e-H/«) . (85) 

Eq. (85) is the Green function of the three-dimensional bi-Helmholtz-bi-Laplace equa- 
tion. It is worth noting that the Green tensor (84) with (85) is in agreement with the 
corresponding expression derived by Zhang et al. (2006). On the other hand, the Green 
function of the bi-Helmholtz equation is given by (e.g., Lazar et al. (2006b)) 

G ( fl » = 4^|)fi( e ^'-^)' (86 > 

In the framework of gradient elasticity of bi-Helmholtz, the differential operator of bi- 
Helmholtz type (77) appears in Eqs. (15)-(18), (23), (27), (30) and (32). 

If we use Eqs. (A. 9) and (A. 10) for the differentiation of Eq. (84), we obtain the explicit 
form of the three-dimensional Green tensor of the bi-Helmholtz-Navier equation 



Gij(R) 



167T/x(l — v) 



2 1 

^cjR 



|( (3 -4, ) (l-- [ l- | ( C? e-^-^e-^)) 

-^U cte ~ R/cl - cte ' R/C2 ) 



I R i R i( 1 6 ( C 1+ C 2) , 6 J_ U e -R/ C1 _ c 4 -R,^ 



r 3 V R 2 ci- ci R 2 

I 6 !. ( J> p-ii/ci _ „3 p ~R/c 2 \ _j_ f ( ,,2 -fl/ci _ 2 p -i?/c 2 

(87) 

The Green tensor (87) gives the non-singular displacement field Ui = Gijfj of the 
Kelvin point force problem, in the framework of gradient elasticity of bi-Helmholtz type. 
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3.1 Dislocation loops 

The calculation of the characteristic fields of a dislocation loop in gradient elasticity of bi- 
Helmholtz type, is analogous to the technique used in gradient elasticity of Helmholtz type. 
The only difference in the results is that now the Green function (86) and the elementary 
function (85) of bi-Helmholtz type enter the characteristic fields of a dislocation loop. In 
gradient elasticity of bi-Helmholtz type, the dislocation density tensor (36) and the plastic 
distortion tensor (38) are given in terms of the Green function of bi-Helmholtz type (86). 
Thus, they are calculated as 



In the limit R — > 0, the integrands of Eqs. (88) and (89) are non-singular at the disloca- 
tion line in contrast to the corresponding ones, Eqs. (37) and (39), calculated in gradient 
elasticity of Helmholtz type. On the other hand, the elastic distortion tensor (45), the 
elastic strain tensor (46), the elastic dilatation (47), the elastic rotation vector (48), the 
stress tensor (49), and the displacement vector (55) are given in terms of the elementary 
function (85) and only (85) has to be substituted in these formulae. The explicit formu- 
lae are not reproduced. The only difference between the fields of a dislocation loop in 
gradient elasticity of bi-Helmholtz type, and of Helmholtz type is that the Green func- 
tion of bi-Helmholtz type (86) and the elementary function (85) have to be substituted 
instead of the Green function of Helmholtz type (28) and the elementary function (17). 
For the derivatives of the function (85), Eqs. (A.7)-(A.12) can be substituted into the 
corresponding formulae. The characteristic fields of a dislocation loop in gradient elas- 
ticity of bi-Helmholtz type retain all the analytical tensor structure of the corresponding 
classical formulae. 

The triple stress tensor of a dislocation loop is easily obtained if the stress tensor is 
substituted into Eq. (71). In gradient elasticity of bi-Helmholtz type the fields produced 
by a dislocation loop are smoother that those predicted by gradient elasticity of Helmholtz 
type. 

3.2 Straight dislocations 

In this subsection, the modified Mura equation (40) is used for gradient elasticity of bi- 
Helmholtz type. The technique of Green functions is used in order to determine the 
non-singular elastic distortion of straight dislocations. 

3.2.1 Screw dislocation 

The Green function of the anti-plane strain problem in gradient elasticity of bi-Helmholtz 
type is the Green function of the two-dimensional bi-Helmholtz-Laplace equation and is 




(88) 



(89) 
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given by (see Eq. (B.34)) 

G ZZ {R) = -^-{lE + In R + — [c\K Q {R/ Cl ) - c 2 2 K (R/c 2 )] } , (90) 

where R = \/(x — x') 2 + (y — y') 2 . The gradient of the Green function (90) is calculated 

as 

G zz , k {R) = - ||{l - ^L^faRK^Rfa) - c.RK^R/c,}}. (91) 

If Eq. (91) and a° zz = b z 5(x)5(y) are substituted into Eq. (40), the elastic distortion 
produced by a screw dislocation with Burgers vector b z is obtained 

Pzx = -^--A l -^-^[cirMr/d) - c 2 rK x {rjc 2 )\ ) , (92) 

■Z7T T v. C-^ C<2 ^ 

P*y = ^- 2 U- -Y^hrKi(r/ Cl ) - c.rK^r/c,)} } , (93) 

where r = ^ x 2 + y 2 . Eqs. (92) and (93) are in agreement with the expressions obtained 
by Lazar and Maugin (2006a). 

The Green function (90) gives the non-singular displacement field u z = —G zz f z of a 
line force with the magnitude f z in the framework of gradient elasticity of bi-Helmholtz 
type. 



3.2.2 Edge dislocation 

The plane strain problem of an edge dislocation is now investigated. The Green tensor 
of the plane strain problem in gradient elasticity of bi-Helmholtz type is found as (see 
Eq. (B.30)) 

Gij(R) = -1T~ SiAlE + \nR + — [cjKo(R/ Cl ) - c 2 2 K (R/c 2 )] } 
+ 16?rM 1 1 _ ;/) 9id 3 {r 2 ( 1e + In R) + 4(c? + 4) ( 7E + In R) 



+ ~2 2 \c\K {R/ Cl ) - 4K (R/c 2 )} } . 

Ci Co L -I * 



(94) 



The gradient of the Green tensor Eq. (94) is calculated as 



Gij,k{R) 
2 



87r/x(l - v) 



fO A \ J R'k r Rj r Ri . n RiRjRk 

(3-4,)^— -^--^- + 2-^5- 



R i R j R k \f2(c 2 + c 2 ) 



~*~ ~p>2 {^ijRk + SikRj + SjkRi ~ 4 ~^ 2 ~ ") ( 



- 4(1 -u)5. 



Rk 



R 2 J \ R 2 
RiRjRk 



r 2 - r 2 
c l c 2 



c\K 2 {R/ Cl ) -c 2 2 K 2 (R/c 2 )^ 



13 R 



R s 



) "1T^> k^l (^M) " C 2^1 
' Ci Co L 



(95) 
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If Eq. (95) and a® z = b x 5(x)5(y) are substituted into Eq. (40), the non- vanishing compo- 
nents of the elastic distortion of an edge dislocation are found as 

- 7273 3V [cirK 1 (r/c 1 )-c 2 rK 1 (r/c 2 )] 

- ^cl-^J K K ^ - ° 2 K ^°^ } ' ^ 
° K X \(3-2u)- 2 4- 4{c2l ] 4) (x 2 -3y 2 ) 



xy ~ 4tt(1 - v) r 2 

2(y 2 + (1 - u)r 2 ) 
(c 2 - c 2 ) r 2 

„2 o„,2^ 



[cirXi(r/ci) - c 2 rK 1 (r/c 2 )] 



+" § [^(r/d) " ^ 2 (r/c 2 )] } , (97) 

I Ci Co ) / J 



(cj - c 2 ) r 

o b * X In o A i V i 4 ( C 1 + C 2) / 2 q 2n 



2(y 2 - (1 - z/)r 2 ) 
+ 772 717772^ [cirK^r/c^ - c 2 rK 1 (r/c 2 )] 



(cj -Co)r 
2{x 2 - 3y 



h LI(l-2v)- — - 4(C ' + ° l) (v 2 - 3x 2 ) 

[cir Ki (r/ ci) - c 2 r (r/ c 2 )] 

+" 73 Iv2 K ^(r/ Cl ) - c\ K 2 (r/c 2 )} } , (99) 

which are in agreement with the formulae given by Lazar and Maugin (2006a). 

The two-dimensional Green function (94) gives the non-singular displacement field, 
Ui = —Gijfj, of a line force with magnitude fj calculated in the framework of gradient 
elasticity of bi-Helmholtz type. 



J yy 



(4- 


c|) r 2 


2v) - 


2x 2 








2{x 2 - 


- vr 2 ) 


r 2 (c\ 


-4) 


2(y 2 - 


-3x 2 ) 


(4- 


c 2 )r 2 



[cjK 2 (r/ Cl )-clK 2 (r/c 2 )}}, (98) 



4 Conclusions 

Non-singular dislocation fields are presented in the framework of gradient elasticity. The 
technique of Green functions is used. The Green tensors of all relevant partial differential 
equations of generalized Navier type were calculated. For the first time, the elastic distor- 
tion, plastic distortion, stress, displacement, and dislocation density of a closed dislocation 
loop, using the theories of gradient elasticity of Helmholtz type and of bi-Helmholtz type 
were calculated. Straight dislocations using Green tensors were revisited. Such gener- 
alized continuum theories allow dislocation core spreading in a straightforward manner. 
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Table 1: Comparison of the basic quantities in different dislocation theories (classical 
dislocation theory, theory of Cai et al. and gradient theory of Helmholtz type) 



Classical theory (DeWit, 1960) 


Cai et al. (2006) 


Lazar [this paper 




R 

AAR = -8n5(x) 


R a 

A A R a = -8nw 
R a = R*w 
w chosen to obtain: 

R a = [ i? 2 + a 2 ] l/2 

a - arbitrary constant 
w = 15a 4 /[87r(r 2 + a 2 ) 7 / 2 ] 


A(R) 
AAA(R) = -8nG 
A(R) =R*G 
LG = 6(x), L = l-£ 2 A 
A(R) = R+?£(l-e- R / e ) 
£ - characteristic length 
G = e- r / e /[An£ 2 r] 



In classical dislocation theory the dislocation function is a Dirac delta function, 5(x), 
without core spreading. In the non-singular approaches by Cai et al. (2006) and Lazar, 
presented in the present paper, the dislocation spreading functions are w and G, respec- 
tively (see Table 1). In the theory of gradient elasticity all formulae are closed in contrast 
to the theory of Cai et al. (2006) where the spreading function w is determined in a so- 
phisticated way in order to obtain R a = [R 2 + a 2 ] 1 ^ 2 . Due to the use of simplified theories 
of gradient elasticity, the dislocation fields retain most of the analytical structure of the 
classical expressions for these quantities but remove the singularity at the dislocation core 
due to the mathematical regularization of the classical singular expressions. In gradient 
elasticity of Helmholtz type, the characteristic length £ takes into account the informa- 
tion from atomistic calculations as discussed in this paper. In a straightforward manner, 
the length £ determines the dislocation core radius. Therefore, in gradient elasticity it 
is not necessary to introduce an artifical core-cutoff radius. It should be mentioned that 
the characteristic lengths which arise in first strain gradient elasticity (e.g., Maranganti 
and Sharma (2007); Shodja and Tehranchi (2010)) and in second strain gradient elastic- 
ity (e.g., Zhang et al. (2006); Shodja et al. (2012)) have been recently computed using 
atomistic approaches. 

The obtained results can be used in computer simulations and numerics of dislocation 
cores, discrete dislocation dynamics, and arbitrary 3D dislocation configurations. The 
results can be implemented in dislocation dynamics codes (finite element implementation, 
technique of fast numerical sums), and compared to atomistic models (e.g., Ghoniem and 
Sun (1999); Li and Wang (2008)). 
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A Appendix: A and its derivatives 



In gradient elasticity theory, the stress tensor, the elastic distortion tensor, the elastic 
strain tensor, and the displacement vector of a dislocation loop are given in terms of 
derivatives of the elementary function A. 



A.l Helmholtz type 

For gradient elasticity of Helmholtz type, the elementary function A is given by 

A - R+2 4 ( 1 - e " s/ ')- 

Higher-order derivatives of A are given by the following set of equations 



Ri 
~R 



1 - 



(l - e"^) 



where R; = X; — x', 



A — y 



RiRj 



11 

+ R e 



6£ 2 



-R/e 



(A.l) 



(A.2) 



l-^(l-e^j + (2 + i je^ 

(A.3) 



A, 



S (l-e' 



-R/e 



). 



(A.4) 



and 



,ijk 



$ij Rk + Sik Rj + Sjk Ri 



R 3 



1 - e- R / £ 



) + ( 2+ i^ 



+ 



3RiRjRk 



R 5 



l--^fl-e-^ 



\ /, 1M 2R\ R/f 

) + ( 4 + ir + 3F) e ' 



A 



iik 



2Rk 



(l-(l + f)e-). 



(A.5) 



(A.6) 



The expressions (A.1)-(A.6) are non-singular. For R — y 0, they are either zero or finite. 



A.2 Bi-Helmholtz type 

In gradient elasticity of bi- Helmholtz type, the elementary function A reads 

2(4 + 4) 2 



A = R + 



R 



cl R 



(A.7) 



The higher-order derivatives of A are given by 



A - Ri 



2(c? + c§) 
R 2 



c\ - c| R? 



Cn e 



cj-c 2 2 R 



2 _ lc 3 e -«/ci_^-il/c 2 



c 2 e 
(A.8) 
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■XI 



R 



RiRj 



1 

+ 



R 2 

6(4 + 4) 

R 2 



c 2 - c 2 R 2 



c(-c 2 R 



+ 



6 



c 2 - c 2 R 2 



— (cW R/ci 



4e~ R/c2 ) + 



C 2 i? 



— ere 



-R/ci 



Ci C 



-R/C! 



Co e 



-r/c 2 \ 



(A.9) 



2 



1 - 



(4 



2 e -R/a _ c 2 e -R/c 2 \ 



(A.10) 



A _ S ij R k + $ik Rj + $jk Rj [ 1 6(4 + 4) i 6 1 fd-R/ Cl n i -R/c 2 \ 



+ 



R 3 
6 J 



c 2 - 4 R 

3RiRjRk 



1 1 11 + 

.^ c 3 e -^i_ c 3 e -i?/ C2 ^ 2 



i? 5 



1 Z 

i? 2 + C 2 -C 2 i? 2 ^ 6 



+ 



R 2 



10 1 ( r ?> p -R/ci _ r 3 p-ii/ca^ , 4 f r 2 p-fl/ci _ r 2 p -ii/c 2 \ 

ci-ci^^i 6 ° 2e ; + C 2_ C 2^ c i e c 2 e ; 

+ 3(^I)( C -^ Cl --^ /C2 )] 



and 



A-,iik 



2Rk , 



* (c?e-^- C ^e-^)- * 



The expressions (A.7)-(A.12) are non-singular. 
(A.12) reduce to Eqs. (A.1)-(A.6). 



(A.12) 

In the limit c 2 ->■ and ci = £, Eqs. (A.7)- 



B Appendix: Green tensors of generalized Navier 
equations 

The following notation is used for the n-dimensional Fourier transform (Guelfand and 
Chilov, 1962) 

/+00 
f(r)e +ikr dr, (B.l) 
-00 

~ 1 r+°° ~ 

f(r) = ^ [/(*)] = 7^ f( k ) ^ rdk ■ (B-2) 
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We have (Wladimirow, 1971; Nowacki, 1986) 



TT-l 

(2) 
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^(3) 


" 1 " 


^(2) 


' 1 " 


^(3) 


" 1 " 


^(2) 


1 




fc 2 + 4 


^(3) 


1 




U 2 + ^J 



— (tb + In v 7 ^ 2 + y 2 ) , 



1 



4vr ^2 + y 2 + z 2 ' 
1 ^^(T^ + lnv^T?), 



8tt 



v^ 2 + y 2 + z<2 , 



exp ( — \/ x 2 + y 2 + z 2 1 c) . 



47Ta/x 2 + y 2 + 

B.l Green tensor of the Helmholtz-Navier equation 

The Green tensor of the Helmholtz-Navier equation is defined by 

(1 - £ 2 A)(M,A + (A + fi)dA)Gij(r) = -Sij6(x) . 

The Fourier transform of Eq. (B.9) reads 

(1 + £ 2 k 2 )(fi6uk 2 + (A + n)kiki)Gij(k) = 5 l3 , 

where A = 2/j,v/(l — 2v) and v is Poisson's ratio. The Fourier transformed 
is found as 



(B.3) 
(B.4) 

(B.5) 
(B.6) 
(B.7) 
(B.8) 
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(B.9) 

(B.10) 
Green tensor 

(B.H) 



Using partial fractions and the inverse Fourier transform, we find 



■y / kjkj 
W\k 4 (l + e 2 k 2 ) 
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-did l; T, 
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(B.12) 



and 



k 2 (l + £ 2 k 2 ) 



1 



,A; 4 (1 + £ 2 k 2 ) 
87T V r r / ' 



(B.13) 
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the three-dimensional Green tensor of the Helmholtz-Navier equation is calculated as 

1 



Gij(r) 



167r//(l — v) 

where r = \J x 2 + y 2 + z 2 . 
On the other hand, using 



2(1 - u)S i:j A - didj 



r + —(\ < 
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B.14) 



F m U 4 (i + £ 2 k 2 ) 
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didjjr w\k\i + e 2 k 2 ' 
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1-2 1,2 1 1 



8tt 



(r 2 ( lE + lnr) + 4£ 2 ( 7£ + lnr) + K (r/£))^j (B.15) 



and 



(^(^(l+W)) ^(k 2 P + i) 



2tt 

the two-dimensional Green tensor of the Helmholtz-Navier equation is obtained as 
1 



(B.16) 



Gijir) = -— 6^ + lnr + K (r/£) } 



^ [r 2 ( lE + In r) + A£ 2 ( lE + In r + K (r /£) ) } , (B. 17) 



167r//(l — 1/) 
where r = ^ x 2 + y 2 . 

B.2 Green function of the Helmholtz-Laplace equation 

For the anti-plane strain problem, the Green tensor of the Navier-Helmholtz equation 
reduces to the Green function of the two-dimensional Helmholtz-Laplace equation which 
is defined by 



(l-£ 2 A)AG zz (r) = — S(x). 



// 



The Fourier transform of Eq. (B.18) reads 

(l + £ 2 k 2 )k 2 G zz (k) 

The Fourier transformed Green function is 

G zz (k) ' ' 



1 

// 



H k 2 (l + £ 2 k 2 ) ' 

Using Eq. (B.16), the two-dimensional Green function is calculated as 

G zz (r) = [ lE + lnr + K (r/£) } . 



(B.18) 



(B.19) 



(B.20) 



(B.21) 
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B.3 Green tensor of the bi-Helmholtz-Navier equation 

The Green tensor of the bi-Helmholtz-Navier equation is defined by 

(1 - c 2 A)(l - 4A)(fi6uA + (A + ^did^djir) = -5, tl 5(x) . 
The Fourier transform of Eq. (B.22) reads 

(1 + c\k 2 )(l + c 2 2 k 2 ){^5 a k 2 + (A + rikik^G^k) = 5i 
The Fourier space Green tensor is 
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A/ M+cf) 
8tt V r 



ct — Co r 



))■ 

(B.26) 



the three-dimensional Green tensor of the bi-Helmholtz-Navier equation is found as 



1 



167T/x(l — Z/) 



2(1 - u)SijA - didj 



2( 



:f - Co r 



4 p-»*/ c i _ c 4 e ^ r / c 2 



(B.27) 



where r = a/x 2 + y 2 + z 2 . 

In two dimensions, we use the formulae 



Jr ( 2 ) 1 ( 



P(l + c 2 A; 2 )(l + c 2 P) 



) T w(k 2 c 2 -c 2 (a; 2 + 4 k 2 + ^)) 

c l c 2 

= "7^(^ + lnr + ^3^2 [c?^o(r/ci) - c^ (r/c 2 )]) 

(B.28) 
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and 

kj k 7 



^ ( kHx + c 2p/ a , C 2 P) ) d A • F (2 ) 1 ( 



fc 4 (l + c?fc 2 ) (1 + clF)/ 1 J ( 2 ) VA; 4 (l + c?P)(l + ciP)- 



3 ( 2 ) Vfc 4 A; 2 c 2_ c 2VA;2 + 4 F + 4 

c l c 2 

r 2 ( 7£ + lnr) + 4(c 2 + c 2 .)^ + lnr) 

[ C ^o(r/ Cl )-c^o(r/ C2 )]). (B.29) 



8tt 

4 



Eventually, the two-dimensional Green tensor of the bi-Helmholtz-Navier equations is 
obtained as 

Gij{r) = -^^^{yE + hir + ^ [c 2 if (r/ci) - 4Ko(t7c 2 )] | 
+ 16^(1 ^{i? 2 (7s + lnr) + 4(c 2 + c 2 ) ( 7b + lnr) 

+ "T^ k^ (r/ Cl ) - c^ (r/c 2 )l } . (B.30) 

B.4 Green function of the bi-Helmholtz-Laplace equation 

For the anti-plane strain problem, the Green tensor of the Navier-Helmholtz equation re- 
duces to the Green function of the two-dimensional bi-Helmholtz-Laplace equation which 
is defined by 

(1 - c 2 A)(l - c 2 2 A)A G zz (r) = -- 5(x) . (B.31) 
The Fourier transform of Eq. (B.31) reads 

(1 + c\k 2 )(l + c\k 2 )k 2 G zz {k) = - . (B.32) 
The Fourier transformed Green function is 

Using Eq. (B.28), the two-dimensional Green function is obtained as 

G zz {r) = ~^{lE + lnr + - J-^ [c 2 K (r/ Cl ) - c 2 2 K (r/c 2 )} } . (B.34) 
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